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Abstract 



We solve a model of interacting electrons coupled to longitudinal phonons 
using the density matrix renormalisation group method. The model is 
parametrised for polyenes. We calculate the ground state, and first excited 
odd-parity singlet and triplet states; and we investigate their energies, and 
bond length changes and fluctuations for up to 30 sites. The transition energy 
and the soliton width of the triplet state show deviations from the adiabatic 
approximation for chain lengths larger than the classical soliton size, because 
of de-pinning by the quantised lattice fluctuations. 
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The inter-play of electron-electron interactions and electron-lattice coupling in polyene 
oligomers and trans-polyacetylene, (CH) X , results in a rich variety of low energy excitations. 
These excitations include triplet states of soliton-antisoliton pairs, singlet states comprising 
bound pairs of triplets, and exciton-polarons. Within the adiabatic (or classical) approx- 
imation [Q] the nature and energy of these excitations is now fairly well understood. A 
parametrised Pariser-Parr-Pople-Peierls model, solved within the adiabatic approximation, 
predicts accurate excitation energies for oligomers of up to 20 or so sites |3[]. However, for 
longer chains there are deviations from the polyacetylene thin film results. These discrepen- 
cies are partly explained by the self-trapping of the excited states by the lattice Q. The 
calculated energies deviate from a linear extrapolation in 1/N as the chain length becomes 
larger than the solitonic structures. We can use this deviation in the energy to estimate an 
upper bound for the self-trapping energy. This deviation is 0.4 eV for the optically allowed 
(l 1 ^") state, 0.7 eV for the even-parity singlet (2 1 Ag) state and 0.3 eV for the lowest lying 
triplet (1 3 B^) Furthermore, a linear extrapolation in 1/N of the short oligomer exper- 
imental values predicts infinite chain energies of the 1 X B~ and 2 1 A+ states close to those 
observed in polyacetylene thin films 0, suggesting that self-trapping may be a partial arte- 
fact of the adiabatic approximation. Thus, the question remains as to the role of quantised 
lattice fluctuations, both on the dimerisation of the ground state, and to the de-pinning of 
the excited states. These fluctuations are the subject of this paper. 

There have been a number of studies of quantised lattice dynamics in the ground state 
of the uncorrelated Su-Schrieffer-Heeger model ||, indicating that fluctuations in the bond 
length are comparable to the bond length changes, but that the Peierls dimerisation is 
stable against such fluctuations. There has also been a variational Monte Carlo study of 
an interacting electron-phonon model J7|. However, there have been no studies of excited 
states, as the incorporation of quantised lattice dynamics into the correlated Pariser-Parr- 
Pople-Peierls model presents a formidable challenge. 

The advent of the density matrix renormalisation group (DMRG) method [JTTJ , (TJj has 



enabled definitive model studies of correlated electron systems, including long range inter- 



actions H and dynamical phonons ||, ||10|| . In this work we report the results of extensive 
calculations on a model system which for the first time affords us insight into the effect of 
quantised lattice dynamics on the properties of excited states of long polyenes. Electrons, 
interacting via long-range Coulomb forces, are coupled to longitudinal phonons. The key 
results of this calculation are that the de-pinning of excited states due to quantum lattice 
fluctuations can become substantial as the conjugation increases. In particular, there is a 
marked reduction in the energy and increase in the soliton width of the triplet excited state. 
The Hamiltonian, with free boundary conditions, is defined as [IU| [H, 

H = huJ2 [ b l b i + -)+hu (blh + tf N b N + l) - hw B t+1 B { + 2Ttg £ {B i+1 - B { ) , 

i=2 ^ ^' i=l i=l 

N / 1\ / 1\ 1 N 

N-l 

- t (1 + 9 {B i+ i - Bi)) (cl +la c ia + clc l+la ) . (1) 

i=l, cr 

b\ (bi) creates (destroys) a phonon and c\ a (c ia ) creates (destroys) an electron on site i. 
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Bi = (&} + bi)/2, g = (\-nhuj /2t) 1 ' 2 and u = V2u = j2K/m. We use the Ohno function 



for the Coulomb interaction: = U/J1 + (f/r^/14.397) 2 , where the bond lengths are in 

A. The single and double bond lengths used in the evaluation of Vij are 1.46A and 1.35A, 
respectively, and the bond angle is 120°. t = 2.539 eV, U = 10.06 eV, T = 0.602, A = 0.115 
and hcu = 0.2 eV |EJ. 

The essential approach we adopt is an extension of the local Hilbert space reduction of 



14| for a representative repeat unit, namely two lattice sites. Once a repeat unit Hilbert 



space is optimised it is then augmented with the system block in the standard finite lattice 
algorithm |TT[. Since the classical lattice geometry of excited states changes as the chain 
length increases, there is no a priori reason to suppose that the optimal repeat unit electron- 
phonon basis for the shortest chain is appropriate for longer chains. Thus, it is generally 
necessary to perform in situ optimisation, i.e. a repeat unit Hilbert space is re-optimised 
when it forms part of the target chain size. Generally, we expect in situ optimisation to be 
necessary whenever the short scale properties are modified by the long scale properties. 

We now outline the procedure in more detail. We begin with a six site lattice, composed 
of three repeat units, and optimise the repeat unit electron-phonon Hilbert space. This is 
done by retaining the optimised states and 'folding-in' some 'bare' electron-phonon states 
(typically 16). Once the full electron-phonon basis has been swept through for a particular 
repeat unit, the same procedure is applied to the next repeat unit until convergence is 
achieved. Next, two repeat units are augmented to form a four site block for the next chain 
size, i.e. 10 sites. The optimised basis for each repeat unit is retained. For 10 sites and 
greater, each repeat unit is re-optimised (including the end units) by sweeping through the 
electron-phonon basis on the first finite lattice sweep. (In situ optimisations in subsequent 
finite lattice sweeps were found to be unnecessary.) During the in situ optimisation only a 
few states (typically < 50) are retained for the environment blocks. After completing the 
sweep the optimised states of the repeat unit are retained for augmentation with the left 
hand block. Typically, 150 states are used for the system and environment blocks during 
augmentation. 

A key goal of this work is to study excited states, which we do by exploiting the particle- 
hole (J) and spin-flip (P) symmetries of Eq. (1) |15|. The inversion symmetry is measured 
at the middle of a finite lattice sweep. We have checked that setting J = +1 and P = +1 
targets the l 1 ^ state, setting J = — 1 and P = +1 targets the 1 X B~ state, and setting 
J = +1 and P = — 1 targets the l 3 B+ state. 

We now turn to the convergence tests. We first establish convergence with respect to 
the number of optimised states per repeat unit. Table I shows the ground state and the 
l l B>~ transition energies for the six site chain for a maximum number of two and three 
bare phonons per site. We see that with 64 states the transition energy has converged to 
within 0.001 eV. Next, we consider the ground state and the 1 X B~ transition energies as a 
function of the maximum number of bare phonons per site, as shown in table II. We see 
that the transition energy has essentially converged to 0.001 eV with five phonons per site, 
and to within 0.02 eV with two phonons (which is better than experimental accuracy). The 
converged l 1 -^ excitation energy of 4.62 eV is very close to the classical result of 4.65 eV. 
Notice also, that the average phonon occupation per site is ca. 0.2. Finally, we consider the 
convergence with super block Hilbert space size for 18 sites. The convergence of the ground 
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state energy is reasonable for up to 180000 states, and the transition energy has converged 
to better than 0.01 eV. 

Next, we consider in what way quantising the lattice degrees of freedom leads to devi- 
ations from the adiabatic approximation. We calculate the classical phonon displacement, 
qi = {huj/ K) 1 ^ 2 {B i ) jl6(| , the bond length distortion and the root mean square fluctuations 
in the bond length. Fig. 1 shows the staggered bond length changes in the ground state 
of an 18 site chain for up to three phonons per site. Also shown is the classical result [F2J. 
In the middle of the chain the phonon calculation is close to the classical result of a bond 
length distortion of ca. 0.05A. However, towards the end of the chain the phonon calculation 
predicts a somewhat larger distortion. We find that in the limit of long chains the relative 
bond length distortion is ca. 0.9, almost independent of the number of phonons per site, 
and close to previous theoretical || and experimental |I7| estimates. We note that the bond 
length fluctuations do not reduce the average bond length change for the linear polyenes 
considered here. This is because the dimerised ground state is not degenerate with respect 
to the state with the bond lengths reversed (in contrast to cyclic polyenes) ||, and thus 
there is no quantum mechanical tunnelling between the two dimerised states. 

Lattice fluctuations do, however, de-pin the self-trapped solitonic structures of the excited 
states. Fig. 2 shows the triplet state transition energy as a function of inverse chain length 
for up to three phonons per site, and for the classical approximation. The energy in the 
classical approximation deviates markedly from a 1/JV behaviour for chain lengths greater 
than 20 sites, as shown in Fig. 2 and especially in [[|. We intrepret this as a result of 
the electronic wavefunction being trapped by the lattice structure. The phonon calculations 
have essentially converged by three phonons per site. In constrast to the classical result, the 
converged phonon calculation shows a much weaker deviation from 1/N behaviour, leading 
to an expected correction of a few tenths of an eV in the infinite chain limit. 

The solitonic structure of the triplet state supports this de-pinning hypothesis. Fig. 3 
(a) and (b) show the soliton structures for 18 and 30 sites, respectively. The position of 
the defect, at roughly the 4th. bond from the center, is roughly the same for both chain 
lengths, and for both the classical and phonon calculations. We expect this, as its position is 
determined by the electronic component of the wavefunction ||. However, we can see that 
for the classical calculation the soliton width is virtually the same for both chain lengths 
(see also Fig. 4 of 0), whereas for the phonon calculation the soliton width is greater in the 
longer chain. This indicates that the coupled electronic and lattice fluctuations lead to an 
increased delocalisation of the wavefunction, and hence to a lower energy. The 2 l A^ state is 
a bound state of two triplets, and it too is self-trapped by the lattice. Although we cannot 
target this state in the current calculation, the above discussion indicates that it will also 
be de-pinned by lattice fluctuations. 

Finally, we consider the optically allowed excitonic (l l B~) state. According to the 
adiabatic approximation [|18|], this state creates a shallow polaronic distortion of the lattice, 
with self-trapping only becoming important for chain lengths longer than ca. 40 sites @. 
Thus, we would not expect lattice fluctuations to play a significant role for the shorter 
chains which we have so far considered. This is confirmed by the excitation energies shown 
in Table IV, indicating that the quantised lattice calculated energies are within 0.1 eV from 
the classical result, and Fig. 4, showing that the quantum and classical polaronic structures 
are virtually identical. 
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In conclusion, an extended DMRG method has been applied to an interacting electron- 
phonon model of polyenes. Quantum lattice fluctuations are shown to play an important role 
in the de-pinning of the self-trapped excited states, leading to corrections to the adiabatic 
approximation, and to an expected reduction of the transition energies of a few tenths of 
an eV for long chains. Thus, a full quantum mechanical treatment of the Pariser-Parr- 
Pople-Peierls model gives remarkably accurate predictions for the excited state energies of 
polyenes. 

M. Yu. L. was supported by the EPSRC (U.K.) (GR/K86343). 
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TABLES 



TABLE I. 


The ground state (l 1 ^) and l 1 ^, 


t transition energies 


(eV) as a function of the 


number of optimised states per repeat unit for the 


six site chain. 




Bare phonons per site Optimised states 






2 


48 


— 99 0539 


4 61 38 


2 


64 


— 99 0578 


4 6087 


2 


80 


—99 0588 


4 6087 


2 


96 


— 29.0596 


4.6085 


2 


144 (Exact) 


-29.0597 


4.6083 


3 


48 


-29.0687 


4.6229 


3 


64 


-29.0756 


4.6184 


3 


80 


-29.0772 


4.6182 


3 


96 


-29.0787 


4.6178 


3 


144 


-29.0793 


4.6173 



TABLE II. The ground state and \^B U transition energies (eV), and the average ground state 
phonon occupation per site for the six site chain. There are 64 optimised states per repeat unit. 


Bare phonons per site 


E^A+) 


E(l 1 B l 




|£?=1<^>1M+ 





-28.6681 




4.353 





1 


-28.9992 




4.577 


0.090 


2 


-29.0578 




4.609 


0.137 


3 


-29.0756 




4.618 


0.164 


4 


-29.0810 




4.622 


0.185 


5 


-29.0832 




4.623 


0.192 


Classical M 






4.646 





TABLE III. The ground state and ^ l B~ transition energies (eV) as a function of the number 
of system block states (m) and superbock Hilbert space size (SBHSS) for the 18 site chain with 
two bare phonons per site. There are 60 optimised states per repeat unit. 



m 


SBHSS 




E0 1 B~)-E(1 1 A+) 


125 


97700 


-89.3167 


3.141 


154 


140248 


-89.3259 


3.145 


173 


184288 


-89.3304 


3.147 
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TABLE IV. The l x B u transition energies (eV) as a function of the number of sites, N. 



AT 




Maximum numb 


er of phonons per site 




Classical [L4j 







1 


2 


3 




10 


3.466 


3.758 


3.837 


3.864 


3.812 


14 


2.984 


3.320 


3.403 


3.424 


3.374 


18 


2.686 


3.061 


3.147 




3.117 


30 


2.240 


2.710 


2.801 




2.786 



s 



FIGURES 



FIG. 1. The ground state staggered bond length change as a function of bond index from the 
center of the chain for the 18 site chain. Circles: one phonon per site, squares: two phonons per site, 
and diamonds: three phonons per site. The dashed line with crosses is the classical approximation. 

FIG. 2. The triplet state (1 3 -B+) transition energy as a function of inverse chain length. The 
same symbols as Fig. 1 (filled diamonds: three phonons per site), with additionally +: no phonons 
per site (i. e. the Pariser-Parr-Pople model). 

FIG. 3. The triplet state staggered bond length change as a function of bond index from the 
center of the chain. (The end bonds are not shown.) The same symbols as Fig. 1. (a) 18 sites and 
(b) 30 sites. 

FIG. 4. The l l B~ state staggered bond length change as a function of bond index from the 
center of the chain for the 30 site chain. The same symbols as Fig. 1. 
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